Astronomy & Astrophysics manuscript no. paper 


© ESO 2012 


February 23, 2012 





Mean-field transport in stratified and/or rotating turbulence 

A. Brandenburg 1,2 , K.-H. Radler 3 , and K. Kernel 1,2 



(N 

o 

(N 

(N 
(N 



6 

(N 
> 

VO 
(N 
(N 

00 

o 



X 



1 NORDITA, AlbaNova University Center, Roslagstullsbacken 23, SE-10691 Stockholm, Sweden 

2 Department of Astronomy, AlbaNova University Center, Stockholm University, SE-10691 Stockholm, Sweden 

3 Astrophysical Institute Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 

February 23, 2012, Revision: 1.168 

ABSTRACT 

Context. The large-scale magnetic fields of stars and galaxies are often described in the framework of mean-field dynamo 
theory. At moderate magnetic Reynolds numbers, the transport coefficients defining the mean electromotive force can 
be determined from simulations. This applies analogously also to passive scalar transport. 

Aims. We investigate the mean electromotive force in the kinematic framework, that is, ignoring the back-reaction of 
the magnetic field on the fluid velocity, under the assumption of axisymmetric turbulence determined by the presence 
of either rotation, density stratification, or both. We use an analogous approach for the mean passive scalar flux. As 
an alternative to convection, we consider forced turbulence in an isothermal layer. When using standard ansatzes, the 
mean magnetic transport is then determined by nine, and the mean passive scalar transport by four coefficients. We 
give results for all these transport coefficients. 

Methods. We use the test-field method and the test-scalar method, where transport coefficients are determined by 
solving sets of equations with properly chosen mean magnetic fields or mean scalars. These methods are adapted to 
mean fields which may depend on all three space coordinates. 

Results. Wc find the anisotropy of turbulent diffusion to be moderate in spite of rapid rotation or strong density 
stratification. Contributions to the mean electromotive force determined by the symmetric part of the gradient tensor 
of the mean magnetic field, which were ignored in several earlier investigations, turn out to be important. In stratified 
rotating turbulence, the a effect is strongly anisotropic, suppre_ssed along the rotation axis on large length scales, but 
strongly enhanced at intermediate length scales. Also the fl x J effect is enhanced at intermediate length scales. The 
turbulent passive scalar diffusivity is typically almost twice as large as the turbulent magnetic diffusivity. Both magnetic 
and passive scalar diffusion are slightly enhanced along the rotation axis, but decreased if there is gravity. 
Conclusions. The test-field and test-scalar methods provide powerful tools for analyzing transport properties of ax- 
isymmetric turbulence. Future applications are proposed ranging from anisotropic turbulence due to the presence of a 
uniform magnetic field to inhomogeneous turbulence where the specific entropy is nonuniform, for example. Some of the 
contributions to the mean electromotive force which have been ignored in several earlier investigations, in particular 
those given by the symmetric part of the gradient tensor of the mean magnetic field, turn out to be of significant 
magnitude. 

Key words, magnetohydrodynamics (MHD) - hydrodynamics - turbulence - Sun: dynamo 



1. Introduction 

Stellar mixing length theory is a rudimentary description of 
turbulent convective energy transport. Th e mixing length 
theory of turbulent transport goes back to Prandtll (|1925l) 



and, in the stellar context, to IVitensel (|l953l ). The simplest 
form of turbulent transport is turbulent diffusion, which 
quantifies the mean flux of a given quantity, e.g., momen- 
tum, concentration of chemicals, specific entropy or mag- 
netic fields, down the gradient of its mean value. In all 
these cases essentially a Fickian diffusion law is established, 
where the turbulent diffusion coefficient is proportional to 
the rms velocity of the turbulent eddies and the effective 
mean free path of the eddies or their correlation length. 

Mean- field theories, which have been elaborated, e.g., 
for the behavior of magnetic fields or of passive scalars in 
turbulent media, go beyond this concept. In the case of 
magnetic fields, the effects of turbulence occur in a mean 
electromotive force, which is related to the mean magnetic 
field and its derivatives in a tensorial fashion. Examples for 
effects described by the mean magnetic field alone, with- 



out s patial derivatives, are the a-effect (jSteenbeck et al 
1966) and t he pumping of mean ma gnetic flux ()Radlei 
19661 H968I: iRoberts fcSowardll 19751): for more infor ma- 
tion on these topics see, e.g., [ <xause fc Radlerl (|l980l ) or 
IBrandenburg fc Subramanianl ( 20051 ). Like wise the mean 
passi ve scalar flux contains a pumping effect (jElperin et all 
1996). In both the magnetic and the passive scalar cases 
turbulent diffusion occurs, which is in general anisotropic. 
The coupling between the mean electromotive force and the 
magnetic field and its derivatives, or mean passive scalar 
flux and the mean scalar and its derivatives, is given by 
turbulent transport coefficients. 

On the analytic level of the theory the determination of 
these transport coefficients is only possible with some ap- 
proximations. The most often used one is the second-order 
correlation approximation (SOCA), which has delivered so 
far many important results. Its applicability is however re- 
stricted to certain ranges of parameters like the magnetic 
Reynolds number or the Peclet number. In spite of this re- 
striction, SOCA is an invaluable tool, because it allows a 
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rigorous treatment within the limits of its applicability. It is 
in particular important for testing numerical methods that 
apply in a wider range. 

In recent years it has become possible to compute the 
full set of turbulent transport coefficients numerically from 
simulations of turbulent flows. The most accurate method 
for t hat is the test-field method (Sch rinner et all [2005, 
120071 ). In addition to the equations describing laminar and 
turbulent flows, one solves a set of evolution equations 
for the small-scale magnetic or scalar fields which result 
from given mean fields, the test fields. By selecting a suffi- 
cient number of independent test fields, one obtains a cor- 
responding number of mean electromotive forces or mean 
scalar fluxes and can then compute in a unique way all the 
associated transport coefficients. 

Most of the applications of the test-field method 
are based on spatial averages that are taken over 
two coordinates. In the magnetic case this approach 
has been applied to a range of differ ent flows i nclud- 
ing isotropic homogeneou s turbulence (|Sur et all [2008; 
IBrandenburg et ail. l2008al). h omogen eous shear flow tur- 
bulence (Brande nburg et all . l2008bl) without and with 



helicity dMitra et all I2009T ). and turbulent convection 
(|Kapylaet all 12009). One of the main results is that in the 
isotropic case, for magnetic Reynolds numbers R m larger 
than unity, the turbulent diffusivity is given by irw^ ms , 
where the correlation time r is, to a good approximation, 
given by t = (u rms fcf ) . Here, u rms is the rms velocity of 
the turbulent small-scale flow and fcf is the wavenumber of 
the energy-carrying eddies. For smaller i? m , the turbulent 
diffusivity grows linearly with R m . Furthermore, if the tur- 
bulence is driven isotropically by polarized waves, the flow 
becomes helical and there is an a effect. In the kinematic 
regime (for weak magnetic fields), the a coefficient is pro- 
portional to uj ■ u, where uj = V x u is the vorticity of the 
small-scale flow, u. In the passive scalar case, test scalars 
are used to determine the transport coefficients. Results 
have been obtained for anisotropic flo ws in the presence 
of ro tation or stro ng magnetic fields ([Brandenburg et ail . 
2009), linear shear dMadarassv fc BrandenburgLl2010[ ). and 
for irrotational flows (Ra dler et all 1201 ID . 

The present paper deals with the magnetic and the pas- 
sive scalar case in the above sense. Its goal is to compute the 
transport coefficients for axisymmetric turbulence, that is, 
turbulence with one preferred direction, given by the pres- 
ence of either rotation or density stratification or, if the 
relevant directions coincide, of both. (Axisymmetric turbu- 
lence can be defined by requiring that any averaged quan- 
tity depending on the turbulent velocity field is invariant 
under any rotation of this field about the preferred axis.) 
Note that a dynamo-generated magnetic field will in gen- 
eral violate the assumption of axisymmetric turbulence. To 
avoid this problem while still being able to investigate the 
general effects arising from only one preferred direction, we 
assume such fields to be weak so as not to affect the assump- 
tion of axisymmetry of the turbulence. An imposed uniform 
magnetic field in the preferred direction would still be al- 
low ed, but this case wi l l not be investigated in this paper; 
see IBrandenburg et al.l (|2009l) for numerical investigations 
of passive scalar transport with a uniform field. 

Except for a few comparison cases, we always con- 
sider flows in a slab between stress-free boundaries. This 
is the simplest example of flows that are non-vanishing on 
the boundary and compatible with axisymmetric turbu- 



lence. To facilitate comparison with earlier work on forced 
turbulence, we consider an isothermal layer even in the 
density-stratified case, i.e., there is no convection, and the 
flow is driven by a prescribed random forcing. This is 
similar to earlier work on f orc ed ho mogeneous turbulence 
([Brandenburg et all l2008alibl . |2009|), but now we will be 
able to address questions regarding vertical pumping as well 
as helicity production and a effect in the presence of rota- 
tion. This setup allows us to isolate effects of density strat- 
ification from those originating from the nonuniformities of 
turbulence intensity and local correlation length. In addi- 
tion to isothermal stratification, we assume an isothermal 
equation of state and thus do not consider an equation for 
the specific entropy. Hence, no Brunt- Valsala oscillations 
can occur. This assumption would need to be relaxed for 
studying turbulent convection, which will be the subject of 
a future investigation. 



2. Mean-field concept in turbulent transport 

2.1. Mean electromotive force 

The evolution of the magnetic field B in an electrically 
conducting fluid is assumed to obey the induction equation, 



— - V x (U x B - V J) . 



(1) 



where U is the velocity and 77 the microscopic magnetic 
diffusivity of the fluid, and J is defined by J = V x B 
(so that J I (Aq with /j,q being the magnetic permeability is 
the electric current density) . We define mean fields as aver- 
ages, assume that the averaging satisfies (exactly or approx- 
imately) the Reynolds rules, and denote averaged quantities 
by overbarsQ The mean magnetic field B is then governed 
by 



— =V x (UxB + £- riJ) 
ot v ' 



(2) 



where £ — u x b is the mean electromotive force resulting 
from the correlation of velocity and magnetic field fluctua- 
tions, u = U — U and b = B — B. 

We focus attention on the mean electromotive force £ 
in cases in which the velocity fluctuations u constitute ax- 
isymmetric turbulence, that is, turbulence with one pre- 
ferred direction, which we describe by the unit vector e. 
Until further notice we accept the traditional assumption 
according to which £ in a given point in space and time 
is a linear homogeneous function of B and its first spatial 
derivatives in this point. Then, £ can be represented in the 
form 



£ = — a±B — («|| — a±)(e ■ B)e — x B 
-P±J - (/3,| - /3 x )(e -J)e-SexJ 
—n±K — (k|| — K±)(e ■ K)e — jie x K 



(3) 



_J_The Reynolds rules imply that F + G = F + G, F = F, 
FG = FG, dF/dx = dF/dx and dF/dt = dF/dt for any fluc- 
tuating quantities F and G. 



2 



A. Brandenburg et al.: Mean-field transport in turbulence 



with nine coefficients a±, an, . . ., /i0 Like J = V x _£?, 

also K is determined by the gradient tensor VB. While 
J is given by its antisymmetric part, K is a vector de- 
fined by K = e ■ (VJ5) S with (V£?) s being the symmetric 
part of VB. A more detailed explanation of © is given 
in Appendix |3 If e is understood as polar vector (for ex- 
ample Vp/|Vp|, where p is the mean mass density), then 
K is axial and 7, /3±, /3n and /i are true scalars, but aj_, 
Q|i, 5, k± and kii pseudoscalars. (Scalars are invariant but 
pseudoscalars change sign if the turbulent velocity field is 
reflected at a point or at a plane containing the preferred 
axis.) Sometimes it is useful to interpret e as an axial vec- 
tor (for example with f2 being an angular velocity). 
Then, K is a polar vector, /3±, /3ii, <5, k±, kh and \i are true 
scalars but ax, «|| and 7 pseudoscalars. 

We may split £ and B into parts £j_ and B± perpen- 
dicular to e and parts £|| and B11 parallel to it. Then Q 
can be written in the form 



-a±B± — 7<> x B± - {3±Jj 
-k±Kj_ — /ie x K± 
-ayS|| -j9yJ|| _ 



5e x J_i_ 



(4) 



K \\ K \\ 



Let us return to In the simple case of homogeneous 
isotropic turbulence we have a± — aii and j3± — (3\\, and 
all remaining coefficients vanish. Then, ^ takes the form 
£ = aB — T] t J with properly defined a and rjt- These two 
coefficients have been determined by test-fie ld calculations 
(|Sur et all l2008t iBrandenburg et all l2008al) . 

In several previous studies of £, more general kinds of 
turbulence (that is, not only axisymmetric turbulence) have 
been considered, but with a less general definition of mean 
fields, which were just horizontal averages. More precisely, 
Cartesian coordinates (x, y, z) were adopted and the aver- 
ages were ta ken over all x and y so tha t they depend on 
z and t only (Br andenburg et ail l2008aHbl ). This definition 
implies remarkable simplifications. Of course, we then have 
J z = 0. Further, there are no non-zero components of V£? 
other than B x z and B y>z , for V • B = requires B Z:Z = 0, 
and these components can be expressed as components of 
J, viz. B XtZ — J y and B VjZ = —J x . (Here and in what 
follows, commas denote partial derivatives.) This again im- 



plies K 



e x J. As a consequence, this definition of 



mean fields reduces ^ to 

-a±B — (aii — a±)(e ■ B)e — je x B 



(5) 



where ft — f3± + hfi and 5^ 



5KX. Of course, a±, 



aii, 7, ft and 5* are independent of x or y. Clearly, (3± 
and /j, as well as S and n± have no longer independent 
meanings. From @ we may conclude that dB z /dt = 0. If 
we restrict ourselves to applications in which B z vanishes 
initially, it does so at all times and the term with aii — a± 
in ([5]) disappears. Then, only the four coefficients a±, 7, ft 
and 8' are of interest. They can be determined by test-field 



2 Note that the signs in front of some individual terms on the 
right-hand side of in particular of those with a± and Qy 
(perpendicular and parallel a effect) as well as 7 (pumping in 
the z direction), may differ from the signs used in other repre- 
sentations. 



calculations using two test fields independent of x and y 
([Brandenburg et all [2008aHbT ) . 

In this paper we go beyond the aforementioned assump- 
tions in the following respects. Firstly, we relax the assump- 
tion that £ in a given point in space is a homogeneous 
function of B and its first spatial derivatives in this point. 
Instead, we admit a non-local connection between £ and B. 
For simplicity, however, we further on assume that £ at a 
given time depends only on B at the same time, thatjs, we 
remain with an instantaneous connection between £ and 
B. This approximation requires that the mean field varies 
slowly on a time scal e much longer than the turnover time 
of the turbulence; see lHubbard fc Brandenburg! (|2009() for a 
more general treatment of rapidly changing fields. Secondly, 
we consider mean fields no longer as averages over all x and 
y. We define B at a point (x, y) in a plane z = const by av- 
eraging over some surroundings of this point in this plane so 
that it still depends on x and y. In that sense we generalize 
© so that 



£{x) = - J{a ± (x,t)B(x-t) 

+ (a||(x,€) - a±(x,$))(e-B(x - £))e 

+ 7 (x,£)exB(x-$) 

+0±(x,£)J(x - £) 

+ (0 {{ (x,t)-0±(x,t))(e-J(x-t))e 

+5(x,£)exJ(x- (6) 

+k±(x,$)k(x-€) 

+ («U (x, £) - K±(x, 0) (e • K(x - £))e 
+n(x,£)exK(x-$)) d 3 £. 

As a consequence of the axisymmetry of the turbulence, 
the coefficients aj_, ay, . . ., fi depend only via £~ + on 
£ x and £ y . We consider them also as symmetric in £ z . The 
integration is over all ^ space. Of course, £, B, J, and K 
may depend on t. For simplicity, however, the argument t 
has been dropped. 

Let us subject © to a Fourier transformation with re- 
spect to £. We define it by 



F(£) = (2tt)- 3 / F(k) exp(ifc • £) d 3 fc . (7) 

Remembering the convolution theorem we obtain 

£{x) = -(2tt)- 3 J (&±(x,kjE(k) 

+ (a|| (x,k) — a±(x, fc)) (e • B(k)j e 
+j{x, fc) e x B(k) 

+h (x, fc) 7(fc) + (/3|, (cc, fc) - /3 ± (x, fc)) (e • 7(h)) e 
+5(x,k)exJ(k) (8) 
+k±(x, k)K(k) + (kh (x, fc) — ii±(x, fc)) (e • if (fc)) e 
+/i(x, fc) e x if(fc) ) cxp(ifc • x) d 3 fc ; 



see IChatteriee et al.l (|201ll) for a corresponding relation in 
the case of horizontally averaged magnetic fields that de- 
pend only on z. Like a±, an, . . ., (j,, the aj_, aii, . . ., /} are 
real quantities. They depend only via k± — (k x + ky) 1 ^ 2 on 
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k x and k y and are symmetric in k z , i.e., depend only via 
fcii = \k z \ on k z . As a±, au, . . ., p are real and symmetric 
in S,x, £,y and £ z we have 



L1U11 ( 

aj_(x,k) = / aj_(x,£) cos k x £ x cos k y t;y cos fc z £ z d 3 £ (9) 

and analogous relations for fiy, ...,Jx. We note that ax, . . ., 
/Et, taken at fc = 0, agree with aj_, . . ., /i in Equation ((3]). 

2.2. Mean passive scalar flux 

There are interesting analogies between turbulent trans- 
port of magnetic flux and that of a passive scalar (cf. 
Radl er et all 1201 ll ). Assume that the evolution of a pas- 
sive scalar C, e.g., the concentration of an admixture in a 
fluid, is given by 



^ = -V • (UC - DVC), 



(10) 



where D is the microscopic (molecular) diffusivity. Then 
the mean scalar C has to satisfy 



Qn 

— = -V -{UC + F-DVC), 



(11) 



where J- — uc is the mean passive scalar flux, u stands 
again for the fluctuations of the velocity and c = C — C 
for the fluctuations of C. Consider again axisymmetric tur- 
bulence with a preferred direction given by the unit vector 
e. Assume that J in a given point in space and time is 
determined by C and its gradient G = VC in this point. 
Then we have 



F = - 7 c C*e - pfG - 



G)e-£°'e x G, (12 



with coefficients 7 C , (3^ and 8 C . If e is a polar vector, 

7 C is a scalar but 8 a pseudoscalar, and if e is an axial 
vector, -f C is a pseudoscalar but 8 C a scalar, while f}±_ and 
are always scalars. We note that V • (8 c e x G) is only 

unequal zero if 8 is not constant but varies in the direction 
of e x G. 

We may split and G into parts J-± and G± perpen- 
dicular to e, and parts !F\\ and Gu parallel to it, and give 
(fH)) the form 

8 c e x Gi 



P±G± 
-f C eC-^G\\ 



(13) 



Let us now relax the assumption that !F in a given point 
in space and time is determined by C and G in this point. 
Analogously to the magnetic case we consider a non-local 
but instantaneous connection between J- and C . Then we 
have 



+ (x, - €)) (e • G(x - €)) e (14) 

+5 c (x,£)exG(x-£)) d 3 £ . 



As ax, • • •) M in the magnetic case, 7 C , /?x\ and (J 
depend only via £~ + on ^ and £ v , and we consider them 



also as symmetric in £ z . The integration is again over all £ 
space. Note that ^F, C, and G may, even if it is not explicitly 
indicated, depend on t. Applying the Fourier transforma- 
tion defined by ([7]) on HH), we arrive at 



-(27T)- 



f(x,fc)eC(fc) 



+^(SB,fc)G(fe) 

+ ( J 8f(x,fc)- i 8f(x,fc))(e.S(fc))e (15) 
+<5 c (a;, fc) e x G(fc) J exp(ifc • a;) d 3 fc , 



where 7^, /3j% /3j? and 8 C are real quantities. They depend 
only via k\ + ky on fc^ and k y , and only via fc|| on k z , and 
they satisfy relations analogous to ©. We note that 7°, 
/^j 7 , /3f?, and 5 C at fc = agree with 7°, Pf> and <5 C 
in (fl 



3. Simulating the turbulence 

We assume that the fluid is compressible and its flow is 
governed by the equations 



BU 

~~Dt 
D/i 

~Dt 



f + g-Vh-2ilxU + p- 1 V • (2vpS) 
-clV -U . 



(16) 



Here, / means a rando m force which prim arily drives 
isotropic turbulence fe.g.. lHaugen et aill2004l ). g the grav- 
itational force, and h the specific enthalpy. An isothermal 
equation of state, p = pc^, has been adopted with a con- 
stant isothermal sound speed c s . In general a fluid flow in 
a rotating system is considered, $7 is the angular velocity 
which defines the Coriolis force. As usual p means the mass 
density, v the kinematic viscosity and S the trace-free rate 
of strain tensor, Sy = \(Ui,j + Uj.i) — h8ij V ■ U. The influ- 
ence of the magnetic field on the fluid motion, that is the 
Lorentz force, is ignored throughout the paper. 

The numerical simulation is carried out in a cubic do- 
main of size L 3 , so the smallest wavenumber is k\ — 2tt/L. 
In most of the cases a density stratification is included with 
g = (0, 0, —g), so the density scale height is H p = c^/g. The 
number of scale heights across the domain is equal to A In p, 
where A denotes the difference of values at the two edges of 
the domain. The forcing is assumed to work with an average 
wavenumber fcf . The scale separation ratio is then given by 
fcf/fci, for which we usually adopt the value 5. This means 
that we have about 5 eddies in each of the three coordinate 
directions. 

The flow inside the considered domain depends on the 
boundary conditions. Unless indicated otherwise we take 
the top and bottom surfaces z = z\ and z = zi with 
Z2 = —z\ = L/2 as stress- free and adopt periodic boundary 
conditions for the other surfaces. 



4. Computing the transport coefficients 

4.1. Test- field method 

In the magnetic case the coefficient s a±, a\\, . . ., p, a re de- 
termined by the test-field method (Schrinne r et al.l . [2005, 
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l2007t iBrandenburg et alll2008ah . This method works with 
a set of test fields B, called B T , and the corresponding 
mean electromotive forces £, called £ T . For the latter we 
have £ T = u x b T , where the 6 T obey 

b T = V x a T 

(>a ' C7xb T + «xB T + (t t xb T )' + ryV 2 a T , (17) 



Of 



with U and u taken from the solutions of p^|) . For the 
boundaries z = const we choose conditions which corre- 
spond to an adjacent perfect conductor, for the x and y 
directions periodic boundary conditions. 
We define four test fields by 

B ls = {B o sxsysz,0,0) , B lc = (B sx sy cz, 0, 0) 

B 2s = (0,0,B sxsysz) , B 2c = (0, 0, B sx sy cz) (18) 

with a constant Bq- Here and in what follows we use the 
abbreviations 



sx = smk x x , cx — cosk x x 
sy = smk y y , cy — cos k y y 
sz = s'mk z z, cz — cos k z z . 



(19) 



We recall that test-fields n eed not to be solenoidal (see 
iSchrinner etail I2005ll200l . 

We denote the mean electromotive forces which corre- 
spond to the test fields CEHJ) by £ ls , £ lc , £ 2s , and £ 2c . With 
the presentation © and relations like © we find 



1 



~Bq (a± sx sy sz — (<5 — — k±)k z sx sy cz) 
-B (7 sx sy sz + (jj± + -fi)k z sx sy cz) 



1 



—p.)k y sxcysz + (8 + — k±)k x cxsysz) 
1 



£ x z = Bo j8y k y sx cy sz 
£l s = -B O {0. 

£ 2 y = B o (0i_ - ^fl)k x cxsysz - (5 + ^k±)k y sxcysz 
£ 2 Z S — — Bo{u\\ sxsysz + k\\k z sxsycz) 

and corresponding relations for £ x c , . . . ,£ 2c , whose right- 
hand sides can be derived from those in ((20]) simply by 
replacing sz and cz by cz and — sz, respectively. 

In view of the assumed axisymmetry of the turbulence, 
we consider a±, an, . . ., \i in what follows as independent 
of x and y but admit a dependence on z. When multiplying 
both sides of the equations (|2"0)) and of the corresponding- 
ones for £ x c , . . . , £ 2c with sx sy, sx cy or cy sy and averag- 
ing over all x and y, we obtain a system of equations, which 
can be solved for a±, an, . . ., fi. The result reads 

& ± = -{b ss (sz£ x s + cz£ x c )) 
a\\ = -(b ss (sz£l s + cz£ 2 z c )) 
7 - -(b ss {sz£ y s +cz£ y c )) 
Px = -i(B ss (cz^ s - sz£ 1 y c ) + B sc (sz£ 2 x s + cz£ 2 x c )) 
= -\{B ss (cz£ l y s - sz£ y c ) - B cs {sz£ 2 v s + cz£ 2 y c )) 
= (B sc (sz£~l s + cz£l c )) (21) 
<5 = \{B ss (cz£ x s - sz£ x c ) - B cs {sz£ 2 x s + cz£ 2 x c )) 





2 \n \czc x 


SZC X ) 


- R sc ( s 7 P 2s _i_ 
D \ azo y 


CZC y )) 


K± = 


-(B ss ( ct£ 1s 


szc x ) 


4- B cs ( S7? 2s 4 


CZC X j) 




-(B ss (cz£ x s 


- ^£ x c ) - 


f 5 sc (sz£ 2s 4 


- CZ£ 2 y C )) 


K|| = 


-(B ss {cz£ 2 z s 


- ^l c )) 






A = 


-(B ss {cz£ y s 


- sz£l c ) 


- B sc {sz£ 2 x s 4 


- cz£ 2 x c )) 


= 


-{B ss {cz£ y s 


- ^£ y c ) - 


f B cs {sz£f^ 


- cz£f)) . 


where 










b ss = 


4 sx sy/Bo , 


B ss = b ss /k z 




B cs = 


4 cx sy/k x Bo 


, B sc = 


Asx cy/kyBo- 





(22) 



The angle brackets indicate averaging over x and y. 
Although the relations (|2ip and (|22|) contain k x , k y and 
k z as independent variables, the a±, an, . . ., jl should vary 

only via kj_ = {k 2 + k 2 ) 1 ^ 2 with k x and k y , and only via ku 
with k z . 



4.2. Test-scalar method 

In the passive-scalar case the coefficients 7 C , /3r/, fi9 , and 
5 C are determined by the test-scalar method with test 
scalars C and the corresponding fluxes T T . For the latter, 
we have !F T — uc T , where c T obeys 



dc T 
at 



-v.( 



Uc 



«C' + (uc T )' - DVc 



(23) 



Again U and u are taken from the solutions of (fT6|) . 

Ts Tc 

We define two test-scalars C and C by 
(20) U — Co sx sy sz , C C — Cq sxsycz, 



(24) 



where Co is a constant and the abbreviations flT9l are used. 
From (jT4"]) we then have 



J-^ = — Co(P±k x cxsysz — 6 c k y sxcysz) 
J- y = —Co0^k y sxcysz + S G k x cxsysz) 
J^ z = — Co(7x sx sysz + /3iFfc z axsy cz) 



(25) 



and analogous relations for J- ' x , . . . ,J r z with sz and cz re- 
placed by cz and —sz, respectively. 

Analogous to the magnetic case, we assume that j c , (3^, 
/3y , and S c are independent of x and y but may depend on 
z. Analogous to ([2"Tj) we find here 

1 C = -{c ss (szT z + czT z )) 

ff[ = -(C cs (szT x + cz^)) = ~(C sc (szT y + czT y )) 

Pf = -(C ss (czT z -szT z )) (26) 

S c = (C sc {szT x + czT x )) = -(C cs {szT y + czTy)) , 

where c ss , C ss , C sc , and C cs are defined like b ss , B ss , B sc , 
and B cs , with Co at the place of Bq- The angle brackets 
indicate again averaging over x and y. Note that j c , 
/Sp and S c should depend only via kj_ = (k 2 + k 2 ) 1 / 2 on 
k x and k y , and only via fey on k z . 
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4.3. Validation using the Roberts flow 

For a validation of our test-field procedure for the deter- 
mination of the coefficients occurring in @ we rely on the 
Roberts flow. We define it here by 

u = uq{— coskyx sin fcoy , sin fcox cosfcoy, 

2/ cos kox cos koy ) , (27) 

with some wavenumber kg and a factor / which charac- 
terizes the ratio of the magnitude of u z to that of u x and 
u y . We further define mean fields as averages over x and 
y with an averaging scale which is much larger than the 
period length 2ir/ko of the flow pattern. When calculating 
the mean electromotive force £ for this flow, we assume 
that it is a linear homogeneous function of B and its first 
spatial derivatives and adopt the second-order correlation 
approximation. Although the Roberts flow is far from being 
axisymmetric, the result for £ can be written in the form 
([3]), and we have 



2rik Q 
ug(l- 



"II ; 
4/ 2 ) 



7 = 



ft_|_ = ft 



= 0, 



i« = 



8 V k 2 ' 
»u(l~4/ 2 ) 



(5 = 



= 2(/3 x -/3||) 



(28) 



It agrees with and ca n be deduced from results reported in 
iRadler et all (|2002allbh . As for the passive scalar case, an 
analogous analytical calculation of the mean scalar flow T 
leads to (TT2"j) with 



0, /9f = 



%Dkl 



2Dkl 



5 C = . 



(29) 



We may proceed from the local connection of £ with B 
and its derivatives considered in ^ to the non-local ones 
given by © or (J5]). As a consequence of the deviation of 
the flow from axisymmetry, we can then no longer justify 
that coefficients like a_i_(0 depend only via + £ 2 on £ x 
and and coefficients like a±(k) only via k± on k x and 
k y . This applies analogously to the connection of T with C 
and its derivatives and to coefficients like and j3±(k). 

A test-field calculation of the coefficients a±, a\\, 

/Et, as well as 7 , . . ., 6 , has been carried out under the 
conditions of the second-order correlation approximation 
with u given by (H7J) and / = l/y/2. Figure Q] shows the 
results obtained for a±, 0±, fiii and /Et, as well as and 
/8y 5 as functions of k±/k{, with fcf = -\/2fco, for two fixed 
ratios k»/k±. In the limit k±/k{ <C 1 these coefficients take 
just the values of a±, j3±, j3\\, /«, Pi and (39 given in ([28]) 
and (j29|) . For larger values of k±/kf, as to be expected, the 
GU_) ^, /3j and /3y depend also on the ratio of k x 

and k„. 
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Fig. 1. The coefficients a±, j3±, $u, and /Et, as well as /?' 
and for the Roberts flow, calculated in the second-order 
correlation approximation, as functions of kj_/k{, where 
kf = V2ko is the effective wavenumber of the flow. Results 
obtained with k x — k y and kn/k± = l/v2 ~ 0.7 or 

fei|/fc_L = l/16-\/2 ss 0.004 are represented by open squares 
and dotted lines or by open diamonds and dashed lines, 
respectively. Results with k x /k y = 0.75 [fcj_ = (3,4, 0)fci] 
or k x /k y — 5 [fej_ = (5, l,0)fei] and = 0.2 are in- 

dicated by open or filled circles, respectively. Orange and 
black symbols correspond to the first and second expres- 
sions for j3± and /x in (f2"Tj) or for ftp in (|26[) . 



4.4. Dimensionless parameters and related issues 

Within the framework of this paper, the coefficients ax, 
an, . . ., fj, as well as a±, am, . . ., jl, and likewise 7 C , /3r_, . . ., 

S c and 7 C , j3j_, . . ., S c , have to be considered as functions 



of several dimensionless parameters. In the magnetic case 
these are the magnetic Reynolds number i? m = u ms /r)k{ 
and the magnetic Prandtl number P m = v /rj, in the pas- 
sive scalar case the Peclet number Pe = u rms /Dkf and 
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the Schmidt number Sc = vjD, further the Mach num- 
ber Ma = u rms /c s , the gravity parameter Gr = g/c^k{, the 
Coriolis number Co = 2f2/u rms fcf, as well as the scale sep- 
aration ratio kf/k\. 

Throughout the rest of the paper we give the coefficients 
a±, an, 7, and j c as well as a±, an, 7, and j c in units of 
Urms/3, the remaining coefficients (3±, . . ., 5 C and . . ., 

5 in units of u Ims /3k{. The numerical calculations deliver 
these coefficients as functions of z and t. To avoid boundary 
effects, we average these results over —2 < k±z < 1 (see 
Figure[3]below). The resulting time series are averaged over 
a range where the results are statistically stationary, i.e., 
there is no trend in the time series. Error bars are defined 
by comparing the maximum departure of an average over 
any one third of the time series with the full time average. 

In the case of isotropic turbulence it has been observed 
that many of the transport coefficients enter an asymp - 
totic regime as soon as R m exceeds unity (|Sur et all [2008) . 
While this should be checked in every new case again (see 
below), it is important to realize that, according to several 
earlier results (see also lBrandenburg et all . [20091) . only val- 
ues of R m below unity are characteristic of the diffusively 
dominated regime, while for i? m exceeding unity the trans- 
port coefficients turn out to be nearly independent of the 
value of R m . 

We are often interested in the limit k±,ku — > 0, in which 
a±, aii, ... S c turn into a±, an ... 6 C . In this limit, how- 
ever, the test fields and test scalars defined by (fT8|) and 
(pM|) vanish. Unless specified otherwise, we approach this 
limit by choosing the smallest possible non-zero \k x \, \k y \ 
and \k z \, that is, by putting k x — k y — k z — k\. 

In the figures of the next section results for a±, an, 

6 are represented. In all cases in which they are considered 
as results for the limit kn — > they are simply denoted 
as a±, an ... S c in the text. 
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Fig. 2. Co dependence of transport coefficients in a model 
with rotation but zero density stratification, R m i=s 9, P m = 
Sc = 1, Gr = 0, fc f /fci = 5. 



5. Results 

5.1. Homogeneous rotating turbulence 

Let us first consider homogeneous turbulence in a rotating 
system, that is, under the influence of the Coriolis force. 
The angular velocity f2 responsible for this force defines 
the preferred direction of the turbulence, e = In this 

case we expect only contributions to the mean electromo- 
tive force £ from a spatially varying mean magnetic field B, 
and contributions to the passive scalar flux !F from a spa- 
tially varying mean passive scalar C. That is, in ([3]) we have 
only the terms with /3j_, /3u, 5, kj_, ku, and /i, and in (|I2p 
only those with /3jJ, @9 , and 5 C . The terms with f3± and (3\\ , 
as well as those with 0^ and (3^ , characterize anisotropic 
mean-field diffus ivities, and tha t with S corresponds to the 
"O x J eff e ct" rtRadlerl. " 
119711 119801 iRadler et al 



"II 



-Co. These coefficients reach maxima at Co w 1. For 



1969^0, H976I: iKrause fc Radlell . 
200.C while the 5 C term van- 



ishes underneath the divergence and is therefore without 
interest. 

Figure [5] shows the dependence of the aforementioned 
coefficients on Co for R m ss Pe ss 9 and kf/ki = 5. The 
values of f3±, j3\\, (3^_ and , which remain finite for Co — > 
0, are always close together. The other four coefficients vary 
linearly with Co as long as Co is small. Specifically, we find 
5 m —0.1 Co, 5 C —Co, as well as k± sa —0.3 Co and 



rapid rotation, |Co| ^ 1, all coefficients approach zero like 
1/Co. In particular, we have (3± w 1.2/Co and the same 
for /3m, Pi, and $9, further k± w -0.5/Co, R\\ « -1.2/Co, 

6 w — 0.3/Co, and 6 « — 0.6/Co. Furthermore, we find 
that, within error bars, a±, an, 7, and j c are indeed zero. 



5.2. Stratified turbulence 

Owing to the presence of boundary conditions at the top 
and bottom of our domain and the lack of scale separa- 
tion for our default choice of kf/k\ = 5, the turbulence is 
in all cases anisotropic, even if gravity is negligible. The 
ratio of the vertical and horizontal velocity components, 
2um/uj_, is no longer, as in the isotropic case, equal to unity. 

For moderate stratification (g/Cgfci f=s 1), not too large \z\, 
and kf/ki = 5, it takes a value of about 0.9. It decreases 
when the ratio k{ / k\ is decreased; see Table [TJ Figure [3] 
shows the z dependence of 2v%/u\. For strong stratifica- 
tion and a high degree of scale separation, e.g. kf/ki = 30, 
the mentioned ratio comes close to unity. Note, however, 
that smaller values of 2u|/u^ can be can be achieved in 
the non-isothermal case when the effects of buoyancy be- 
come important. 
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Fig. 3. Anisotropy 2u|/it^_ of nonrotating turbulence for 

different stratifications, g/c*ki, and different degrees of 
scale separation, fcf/fci. 

5.2.1. Stratified nonrotating turbulence 

For axisymmetric turbulence in a nonrotating system show- 
ing any kind of stratification in the representation ([3]) of £ 
only the four coefficients 7, /3±, /3n, and p can be non-zero. 
Likewise, in the representation (|12p of J- only the three co- 
efficients 7 C , /3j% and can be non-zero. Figure 2] shows 
their dependence on Gr. It appears that 7 is always close 
to zero, while j c shows a linear increase for not too strong 
gravity. At the same time, (3±, /?||, /3±, and fi9 remain ap- 
proximately constant. We find that p is negative and its 
modulus is mildly increasing with increasing stratification, 
but the error bars are large. 

5.2.2. Stratified rotating turbulence 

For turbulence under the influence of gravity and rotation, 
all nine coefficients a±, . . ., p are in general non-zero, as 
well as all four coefficients j c , . . ., S c . If both gravity and 
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Fig. 4. Gr dependence of the transport coefficients in a 
model with density stratification but zero rotation, P m = 
Sc = 1, R m « 22, Co = 0, fcf/fci = 5. 



rotation are so small that 5 is linear in g and f2, more 
precisely £ contains <7 m f2™, where n and m mean integers, 
only with n + m < 1, aj_ and ay vanish but 7, /3j_, 5 and 
k_i_ may well be unequal to zero. If n + m < 2, all nine 
coefficients may indeed be non-zero. 

Results for stratified rotating turbulence are shown in 
Figure [SJ The error bars are now bigger than either with 
just rotation or just stratification. For Co —> 0, the coef- 
ficients /3±, /3||, p, /3±, (3^ and Sc remain finite. As Co is 
increased, their moduli show some decline. On the other 
hand the moduli of a±, an, 7, 5, k±, km and 7 C increase 
with Co as long as it is smaller than some value below unity 
but decrease again for larger Co. Both a± and otu are neg- 
ative, which is expected for g and Jl being antiparallel to 
each other. Interestingly, /i is finite for small values of Co, 
in agreement with the result when there is only stratifica- 
tion (Figure S]), but with a modest amount of rotation, /1 
is suppressed and grows only when Co has reached values 
around unity. 



5.3. Wavenumber dependence 



in the limit fc 



Ifcl 



0, that is, k±,k\\ —> 0. However, 



Table 1. Dependence of the density contrast Phot /pto P and So far we have cons i der ed the coefficients &±, fij 
the degree of anisotropy 2n|/ttj_, for three different values 
of fcf/fci, on the density stratification g/c^ki for nonrotating 
turbulence. The values of 2m|/u^_ have been obtained as 
averages over the range — 2 < k\Z < 1. 



g/c%ki 


Pbot/Ptop 


fc f = 1.5fci 


1u\ju\ 
fcf = 5fci 


fcf = 30fci 








0.84 


0.99 


1.00 


0.5 


23 


0.84 


0.97 


1.00 


1 


540 


0.66 


0.90 


0.99 



their behavior for larger k, in particular for k up to sev- 
eral fcf, is of interest, too. Most of them decrease like k~ 2 
as fc grows and can be fitted to a Lorentzian profile, as 
has been fo und in earlier calculation using the test-field 
method; see iBrandenburg et al.l (|2008afi . where in fact the 
dependence on fcii was considered. Even earlier work that 
was n ot based on the test-field method showed a declinin g 
trend (|Miesch et all . l200fA IBrandenburg fc Sokoloffl l2002h . 
Nevertheless, as is shown in Figure [6j there are also some 
coefficients that first increase with fc|| , have a maximum 
near fcii = fcf and only then decrease with growing fc|| . 
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Fig. 5. Co dependence of transport coefficients in a model 
with rotation and density stratification, P m = Sc = 1, 
R m w fO, Gr « 0.16, k { /h = 5. 



Fig. 6. fcii dependence of transport coefficients in a model 
with rotation and density stratification, kj_ = v2fci, P m = 
Sc = 1, R m = 12, Co = 1.0, Gr = 0.16, kf/fo = 5. 



Examples for such a behavior are a\\, 5, and k±, while Kti 
peaks slightly below ku — 0.5kf. 

The dependence of the coefficients under discussion on 
k± is shown in Figure [7] Note that our test fields vanish for 
k± = 0, so no values are shown for this case. Note also that 
—&\\, arL d — /tii, which have maxima for ku/ki rs 1 or 
fcii/fef f=a 0.5, show a clear monotonic decline with k±. Only 
— kj_ has maxima with respect to fc||/fcf and fcj_/fcf. 

Most of the results presented in Figure [7] have been cal- 
culated with k x = k y , a. few single ones for aj_ , /3j_ , k± and 

also with k x jky = 0.75 and k x /k y — 0.2. While the re- 
sults for j3± and agree well for all these values of k x /k y , 
there are significant discrepancies with a± and fix- 

5.4. Dependencies on R m and Pc 

Let us finally consider the dependence of all transport coef- 
ficients on i? m or Pe for a case where they are all expected 
to be finite. Therefore we choose again the case with Co = 1 
and Gr = 0.16, which was also considered in Figures [5j-[71 
and keep P m — Sc = 1. 

The results are shown in Figure [5] As expected, some 
of the quantities increase approximately linearly with R m 



if R m < 1, or with Pe if Pe < 1, and seem to level off to 
constant values for larger values of R m , or P m , although 
the uncertainty tends to increase significantly. 

6. Conclusions 

In this paper we have dealt with the mean electromotive 
force and the mean passive scalar flux in axisymmctric tur- 
bulence and have calculated the transport coefficients that 
define these quantities. Unlike most of the earlier work, we 
have no longer assumed that mean fields are defined as pla- 
nar averages but admit a dependence on all three space 
coordinates. The number of test fields and test scalars is 
the same (4 and 2, respectively) as in earlier work using 
planar averages, so the computational cost is unchanged. 

We may conclude from general symmetry considerations 
that the mean electromotive force £ has altogether nine 
contributions: three defined by the mean magnetic field B, 
three by the mean current density J, and three by the vec- 
tor K, which is the projection of the symmetric part of the 
gradient tensor V£? of the magnetic_field on the preferred 
direction. In many representations of £ the last three contri- 
butions have been ignored. Our results underline that this 
simplification is in general not justified. The corresponding 
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Fig. 7. Same as Figure [51 but k± dependence, fcii = k\. 
The filled and open circles denote results for a±, pj_, k±, 
and obtained with k x /k y = 0.75 [k± — (3,4,0)fci] and 
k x /ky = 0.2 [fcj_ = (f,5,0)fci], respectively. 



coefficients k|| and /i are in general not small compared 
to /3j_, and 5. 

It has been known since long that a stratification of 
the turbulence intensity, th at is, a gradie nt of u 2 , caus es a 
pumping of magnetic flux (|Radlerfll966l 119681 Il969bl) . It 
remained however uncertain whether the same effect occurs 
if a preferred direction is given by a gradient of the mean 
mass density ~p while the turbulence intensity is spatially 
constant. In our calculations, which correspond to this as- 
sumption, the value of 7 is not clearly different from zero. 
This suggests that a gradient of the mass density alone 
is not sufficie nt for pumping, what is also in agreement 
with results of lBrandenburg et al.l (|201ll) . This is even more 
remarkable as the corresponding coefficient 7 C which de- 
scribes the transport of a mean passive scalar is noticeably 
different from zer o. Pumping down th e density gradient is 
indeed expected (|Elperin et al.l . |1995|) . An explanation of 
these results would be very desirable. 

In homogeneous rotating turbulence, apart from an 
anisotr opy of the mean -field conductivity, the f2 x J effect 
occurs (Radlcr, I969a ; b). In the passive scalar case again 
an anisotropy of the mean diffusivity is possible. Even if the 



Fig. 8. Dependencies of the transport coefficients on R m 
or Pe in a model with rotation and density stratification, 
P m = Sc = I, Co = 1.0, Gr = 0.16, h/fa = 5. 



flux proportional to l~i x VC is non-zero, it cannot influence 
C, ' 

Let us turn to the induction effects described by K. If 
the preferred direction is given by a polar vector, the corre- 
sponding contribution to the mean electromotive force can 
only be proportional to e x K. We found such a contribution 
in the case of the Roberts flow and also, for turbulence sub- 
ject the Coriolis force, in the results presented in Figure [5] 
and Figures 0HZ1 

Contributions to the mean electromotive force as de- 
scribed here by K occur also in earlier calculations, 



e.g. iKitchatinov et al.l (|l994h or iRtidiger fc Brandenburg! 
(|1995E As a consequence of other notations, however, this is 
not al ways obvious. For example, IRtidiger fc Brandenburg] 
(|l995l) consider a mean electromotive force of the form 



£ = —r)\\J + (rj\\ - Vt)(zJ z - z X VB Z ) 



(30) 



with two coefficients rj^ and 77T (equation (18) of their pa- 
per with ^oJ^'m the sense of the definition introduced here, 
replaced by J; z is our e). It is equivalent to our represen- 
tations Q or (01 of £ if we put there j3± = ^(^ + 7/t), 
/3|l = r/T, M = ~ Vt and all other coefficients equal to 
zero. This implies j3±_ — f3\\ = /i/2, which is in agreement 
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with the relation for fi in equation ([25)1 for the Roberts 
flow. The latter equality is also approximately obeyed for 
turbulence in the presence of rotation, stratification, and 
both; see Figs.^JHl and[Sl respectively. 

If there is moderate rotation (Co s» 1), but no strati- 
fication, we have j3± > /3n; see Figure [2j This means, e.g., 
that for a magnetic field without a component in the di- 
rection of the rotation axis the diffusion along this axis is 
enhanced compared with that in the perpendicular direc- 
tion. In the passive scalar case we have (3^ > , which im- 
plies that the diffusion along the rotation axis is enhanced, 
too. However, stratification enlarges /3ii —j3± and diminishes 

~ P± so t na ^ the diffusion along the rotational axis is 
decreased in both cases considered. In the presence of ro- 
tation and density stratification all three contributions to 
the mean electromotive force described by K are in gen- 
eral non-zero. Here, is smaller than There is now 
also an a effect, which is necessarily anisotropic, and |aJ 
is typically only half as big as |a_i_|; see Figure [5] 

The present work is applicable to investigations of stel- 
lar convection either with or without rotation, and it would 
provide a more comprehensive description of turbulent 
transport properties than what has been available so far 
(Kap yla, et all |2009[ ). The methods utilized in this paper 
can be extended to a large class of phenomena in which 
turbulence with just one preferred direction plays an im- 
portant role. Examples for that include turbulence under 
the influence of a strong magnetic field and/or an exter- 
nally applied electric field leading to a c urren t perm eating 
the system. Turbulence generated by the iBelll (|2004|) insta- 
bility is an example. In addition to density stratification, 
there can be a systematic variation of the turbulence in- 
tensity in one direction. A further example is entropy in- 
homogeneity combined with gravity giving rise to Brunt- 
Valsala oscillations. Pumping effects als o exist in homoge- 
neous flows if the turbule nce is helical (|Mitra et all . l2009t 
iRogachevskii et ail 1201 ll) . By contrast, shear problems or 
other types of problems with two or more preferred direc- 
tions that are inclined to each other (e.g., turbulence in a 
local domain of a rotating stratified shell at latitudes differ- 
ent from the two poles) are not amenable to such a study. 
Of course, although we refer here to axisymmetric turbu- 
lence, problems in axisymmetric cylindrical geometry are 
also not amenable to this method, because the turbulence 
must be homogeneous in one plane. 
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Appendix A: Derivation of relation ([3]) 

We start from the aforementioned assumption according to 
which £ is linear and homogeneous in B and its first spatial 
derivatives, 

Si = aijBj + b m (VB) jk . ( A. 1 ) 



Here aij and bijk are tensors determined by the fluid flow. 
The gradient tensor (VB)ji can be split into an antisym- 
metric part, which can be expressed by J, and a symmetric 
part (VB)^. Therefore we may also write 

Si = n,,7l, - bijJj ~ c ljk (VB)% (A.2) 

with new tensors bij and Cij k , the latter being symmetric 
in j and k. From the further assumption that the flow con- 
stitutes an axisymmetric turbulence we may conclude that 
dij, bij and Cijjk are axisymmetric tensors. Defining the 
preferred direction by the unit vector e we then have 

dij = ai5i.j + a 2 (-ijiei + aze^j , 
b^ = biSi-j + b 2 ei-jiii + b-ihe-j , 
Cijk = ciSjkii + c 2 (5ije k + S ik ej) (A.3) 
+C3(eijieie k + tikihzj) + c A e.ieje k , 

with coefficients ax, &2, • ■ ■, c& determined by the fluid flow. 
Taking (|A.2[) and (|A.3[) together and considering that 

(S tj e k + S tk ej)(VB)% = 2K t , 
(eijieie k + eikieiej)(VB)% = -2(e x K) t , (A.4) 
e.ieje.k{V B) s jk = (e • K)ii , 

we find 

£ = a\B — a 2 e x B — a^(e ■ B)e 

+bxJ - b 2 e xj - b 3 (e -J)e (A. 5) 

+2c 2 ~K - 2c 3 e x ~K + c 4 (e • ~K)e . 

Since (VB),, = there is no contribution with cx- With a 
proper renaming of the coefficients (IA.5I) turns into ([3]). 
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